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We introduce a novel optimization method to design the first examples of photonic quasicrystals 
with substantial, complete photonic band gaps (PBGs): that is, a range of frequencies over which 
electromagnetic wave propagation is forbidden for all directions and polarizations. The method 
can be applied to photonic quasicrystals with arbitrary rotational symmetry; here, we illustrate the 
results for 5- and 8-fold symmetric quasicrystals. The optimized band gaps are highly isotropic, 
which may offer advantages over photonic crystals for certain applications. 
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I. INTRODUCTION 

Photonic materials with complete photonic band gaps 
(PBGs), i.e., frequency ranges over which electromag- 
netic wave propagation is prohibited for all directions 
and polarizations, are artificial dielectric heterostructures 
that enable control of the generation and flow of light. 
Their development 1,2 has evolved dramatically in the 
last decade and their unusual properties have found a 
wide range of applications including efficient radiation 
sources^ telecommunications devices £ sensor s^ and op- 
tical computer chips & Until recently, the only known ma- 
terials with complete PBGs were photonic crystals com- 
prised of a periodic arrangements of dielectric materials. 
In photonic crystals, the opening of the PBG is under- 
stood to be governed by the synergetic interplay between 
Bragg scattering resonances of the periodic dielectric ar- 
ray and the Mie resonances of individual dielectric scat- 
tering centers^ 

A similar combination of scattering mechanisms is 
possible in photonic quasicrystals, in which the dielec- 
tric materials are arranged in a pattern with long-range 
quasiperiodic translational order and rotational sym- 
metries forbidden to crystals (such as 5-fold symmetry 
in two dimensions and icosahedral symmetry in three 
dimensions)^ Quasiperiodic order also produces Bragg 
scattering^. In fact, there have been numerous studies of 
photonic quasicrystals with PBGs in the literature *2r— 
Their band gaps are found to be considerably more 
isotropic and, for two-dimensional quasicrystals, the di- 
electric contrast required to open the band gaps for TM 
polarization (electric field oscillating out of the plane) 
is smaller than the contrast required for the periodic 
counterparts^ Most of the photonic quasicrystalline 
structures considered only have band gaps for either TM 
or TE polarized radiation ^ 13 ! 15 Two examples of com- 
plete band gaps have been found previously ; 12 ^ 6 here 
we present a systematic method to produce substantially 
greater fundamental complete band gaps for photonic 
quasicrystals with arbitrary symmetry. 

In general, the architectures required for optimal TM 
band gaps are quite different than optimal architectures 




FIG. 1: Protocol described in Sec. Ullmaps the vertices of a 
rhombus tiling (red/grey points) into a network of cells whose 
vertices are trivalent (black segments and points). To con- 
struct a photonic material with a complete band gap, edges 
are replaced by a wall of dielectric of finite thickness and the 
vertices are replaced by cylinders of finite radius. 



for TE polarization (electric field oscillating in the plane). 
Structures based on a distribution of isolated dielectric 
inclusions (e.g., rods) is best for TM band baps, while 
structures based on connected dielectric networks are op- 
timal for TE*i£ For the case of quasicrystalline systems, 
with their rich variety of local environments, the problem 
of finding an optimal compromise becomes more difficult. 
Consequently, until now, there has not been a systematic 
procedure for designing quasiperiodic photonic structures 
with sizable PBGs for all directions of propagation and 
polarizations. 

In this paper, we introduce a novel optimization 
scheme discussed in Ref. 18 (where it was applied to 
disordered heterostructures) to design two-dimensional 
photonic quasicrystals with substantial complete PBGs, 
comparable to the largest band gaps found for photonic 
crystals. The examples illustrated here are based on the 
vertices of a Penrose tiling^ a 5-fold symmetric pattern 








FIG. 2: Optimal photonic crystal and quasicrystal heterostructures and their diffraction patterns derived from: (a) a periodic 
6-fold symmetric point pattern; (b) a 5-fold symmetric quasicrystal (Penrose) tiling; (c) an 8-fold symmetric quasicrystal tiling. 



composed of obtuse and acute rhombi and on the ver- 
tices of an octagonal tiling^ an 8-fold symmetric pattern 
composed of squares and rhombi. 

In Section [TTJ we describe how the quasiperiodic pat- 
terns are constructed. In Section IIII[ we present the re- 
sults of the band structure calculations and analyze the 
properties of the resulting PBGs including the field dis- 
tribution of the photonic modes that define the photonic 
band edges. Section llVl provides concluding remarks. 



II. NEARLY OPTIMAL PHOTONIC BAND BAP 

STRUCTURES FROM QUASIPERIODIC POINT 

PATTERNS 



The Penrose and octagonal quasicrystalline point pat- 
terns considered in this work are constructed by the pro- 
jection and cut method^ which consists of projecting 
the points of an n-dimensional hypercubic lattice into 
the three dimensional space, where n = 4 for the octag- 
onal lattice and n = 5 for the Penrose lattice. In gen- 
eral, the lack of spatial periodicity makes it impossible to 
define a unit cell for band gap calculations; we address 
this by using a series of rational periodic approximants 
of the quasicrystal. In the case of the Penrose lattice 
periodic approximants are constructed by the same pro- 
jection method by replacing by a Fibonacci ratio (r n = 
Pn/Qn — 1/1,2/1,3/2,5/3,...). The unit cell is rectan- 
gularly shaped and its area grows as n increases and the 
rational approximant approaches r. In particular, for a 
rational approximant r n = p n /q ni the unit cell has di- 



mensions L x = h(p n + q n (r - l),L y = V3 - r(p n r + q n ) 
and contains N = 10p n (p + 2q n ) vertices of the rhombic 
pattern. In the case of the octagonal quasilattices, we 
replace by its continuous fraction series (p n = p n /q n = 
1/1, 3/2, 7/5, ...), and the resulting periodic approximant 
has a square shaped unit cell of side L = p n + pq n con- 
taining N = (p n + (p 2 + p)q n )(Pn + (p 2 - p)q n ) vertices 
of the square-rhombic pattern. The periodic approxi- 
mants constructed this way results in tilings consisting 
of identical tiles as the original quasiperiodic tiling and 
transition from quasiperiodicity to periodicity is accom- 
modated through the flipping of a certain number of 
tiles. These rational periodic approximants are known 
to be the best approximants of the quasiperiodic lattice 
as they achieve the selected periodicity by introducing 
the minimum density of defects with respect to the per- 
fect quasiperiodic tiling^ The vertices of the tilings ob- 
tained by this method form a sequence of quasicrystal 
approximant point pattern. 

The next step is to find the arrangement of dielectric 
material around the point pattern that produces the op- 
timal TM, TE and complete PBGs. Identifying these op- 
timal dielectric distributions is well known to be a daunt- 
ing computational task, despite the recent development 
of rigorous mathematical optimization methods^"— 

For the case of TM radiation only, a nearly optimal 
decoration for a given point pattern is obtained by plac- 
ing identical dielectric cylinders centered at each point 
and adjusting the radius. For the optimal radius, there 
is a coincidence of Mie and Bragg scattering effects that 
lead to substantial band gaps^ The dielectric cylinders 
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FIG. 3: The fractional PBG width for (a) 5-fold symmetric 
and (b) 8-fold symmetric photonic quasicrystal as a function 
of the number of particles in the periodic approximant for the 
optimal TM (red circles), TE (orange squares) and complete 
(green diamonds) band gaps. 



support Mie scattering resonances and, for frequencies 
above the lowest order resonance, the scattered radiation 
is out of phase with the incident wave. This destructive 
interference prohibits radiation propagation and favors 
the opening of a complete TM gap^ At the same time, 
the quasiperiodic long range order in quasicrystalline sys- 
tems results in a dense (no minimum separation) collec- 
tion of Bragg peaks; however, it is notable that most of 
Bragg peak intensities are infinitesimally small. The first 
set of intense Bragg peaks is associated with scattering 
on planes separated by the average distance in the qua- 
sicrystalline point pattern and can be employed to define 
an effective Brillouin zone. Therefore, as for photonic 
crystals, Bragg scattering also contributes to the PBG 
formation in photonic quasicrystals such that, whenever 
the wavevector of the incident radiation is directed along 



the effective Brillouin zone boundary, the reflected and 
refracted waves interfere to cancel the incoming wave and 
prevent its propagation inside the structure. An optimal 
band gap occurs when the two scattering mechanism re- 
inforce each other on a given spectral range. This ap- 
proach can be applied to find near optimal TM band gap 
for the 5- and 8-fold quasicrystal patterns obtained by 
projection^ 

To obtain the optimal PBGs for TE radiation, where 
the electric field is oriented in the plane of the scatters, 
necessitates a different dielectric arrangement. Instead 
of isolated cylinders, a connected network of dielectric 
with air pockets in between is favored. For example, 
a commonly used configuration for photonic crystals is 
"inverted" compared to the optimal structures for TM 
PBGs, i.e., placing an identical air cylinder at each point 
so as to produce a connected network of dielectric mate- 
rial. However, in the case of quasicrystalline structures, 
we find this method fails to produce sizable fundamental 
TE PBGs. The main reason is that the inverted structure 
has a very non-uniform distribution of dielectric scatter- 
ing regions that broadens the distribution of resonances. 
Similarly, placing walls along the edges of each Penrose 
tile as in Ref. [9j has rather non-uniform connectivity and 
fails to produce an optimal TE PBG. 

To overcome this problem, we employ a protocol in- 
troduced in Ref. [l8| that converts a general point pattern 
into photonic heterostructures with sizeable, nearly op- 
timal TE PBGs. We briefly review the protocol here: 
Consider an arbitrary point pattern denoted by red disks 
in Fig. [TJ In the example in the figure, the point pattern 
is enclosed in a square domain with periodic boundary 
conditions. First, we perform a Delaunay tessellation 
of the original point pattern. This provides a triangu- 
lar partitioning that minimizes the standard deviation of 
the triangle angles around the mean of 7r/6. Next, asso- 
ciated with each triangle is a centroidal point; these are 
connected with line segments to form cells around the 
original points. The nearly optimal heterostructure for 
TE band gaps is obtained by decorating the edges of the 
cell network with walls (along the azimuthal direction) 
of dielectric material of uniform width w, as displayed 
in Fig. [TJ We note that this procedure is universal and 
can be used to generate connected network architectures 
based on periodic, quasiperiodic and disordered point 
patterns. Finally, the optimal structure is obtained by 
varying a single parameter, the cell thickness, and iden- 
tifying the value that maximizes the TE band gap. 

The total number of cells in the network is the same 
as the number of points in the original point pattern and 
each vertex of the network is trivalent. The cells can be 
regarded as individual scattering objects supporting elec- 
tromagnetic resonances, which become coupled to each 
other when placed into the connected network. 

With an additional step, the protocol can be used to 
a heterostructure with a nearly optimal complete PBG, 
i.e., maximal overlap between TM and TE gaps. Namely, 
each vertex of the trihedral network constructed above is 
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FIG. 4: Photonic band structure and density of states for 
5-fold and 8-fold symmetric photonic quasicrystals. (a) Pho- 
tonic band structure and the corresponding density of states 
for a network constructed via the protocol in Sec. [XT] from a 
2/1 periodic approximant of the Penrose tiling with optimized 
cylinder radius r /a = 0.157 and wall thickness w/a — 0.042. 
The structure has a complete PBG of 15.14% with central 
frequency ouc/(2ira) — 0.43. (b) Photonic band structure and 
the corresponding density of states for a network constructed 
from a 3/2 periodic approximant of the octagonal quasilattice 
with optimized cylinder radius r/a=0.1613 and wall thickness 
w/a = 0.0136. The structure has a complete PBG of 13% 
with central frequency ujc/(2ixa) — 0.44. 



decorated with a cylinder or radius r. The number of 
cylinder is equal to twice the number of cells. The opti- 
mal complete band gap obtained by varying the two free 
parameters, the cell wall thickness w and cylinder radius 
r. As shown in Fig. [2j the Fourier spectrum of the deco- 
rated photonic structures obtained through this protocol 
fully preserve the symmetries of the original point pat- 
tern. In comparison, alternative designs based on "in- 
verted" architectures or decorations of the quasicrystal 
tile edges with dielectric walls and its vertices with di- 
electric cylinders fails to open sizable complete PBGs. 



III. PHOTONIC BAND GAP RESULTS 

We compute the photonic band structures for the opti- 
mal heterostructures (described in detail below) using a 
conventional plane- wave expansion formalism^ For this 
approach, which assumes a periodic structure, we use a 
sequence of increasingly accurate periodic approximants 
and check for convergence. For the purposes of illustra- 



tion, we assume the photonic materials are composed of 
silicon (with dielectric constant e = 11.56) and air. In all 
the numerical simulations, we use 32 2 x Np plane waves 
to achieve convergence accuracy of better than 1% for 
the lowest Np photonic bands. Here Np represents the 
number of scattering centers (TM case), cells (TE case) 
or centers and cells (TM+TE case). Most of the band 
structure calculations are preformed around a contour 
along the first Brillouin zone of the respective system, 
which includes the high symmetry k-space points, T = 0, 
X = bi/2, M = (bi + b 2 )/2, and Y = b 2 /2, where 
bi and b 2 are the basis vectors of the reciprocal lattice 
considered. 

TM Band Gaps: To obtain optimal TM PBGs for a given 
symmetry, dielectric cylinders are placed at the vertices 
of the tilings obtained by projection. 

For example, an optimized 5-fold symmetric quasicrys- 
talline pattern derived using the protocol from a n = 5/3 
periodic approximant of the Penrose pattern displays a 
TM PBG of Auj/uc = 39%. An analogous optimized 
structure based on the n = 7/5 periodic approximant of 
the octagonal lattice displays a TM PBG of Auj/ujc = 
42%, where uoc is the gap central frequency. The optimal 
radius of the cylindrical inclusion is r/a = 0.177 for the 
Penrose tiling and r/a = 0.189 for the octagonal lattice, 
were a is the side length of the rhombic tiles obtained by 
projection. As a check, we perform a convergence analy- 
sis as the order of the periodic approximant increases; as 
shown in Fig. [3l the results are fully convergent. 

We find that the band gap for the optimal structures 
always occurs between bands Np and Np + 1 where Np is 
the number of points in the periodic approximant point 
pattern. This is an analogous behavior to the band fold- 
ing that occurs in periodic structures when analyzed with 
a super-cell approach. 

TE Band Gaps: The optimal TE band gap designs are 
obtained by the protocol described in SecJIIl They corre- 
spond to the network of walls connecting the centroids of 
the Delaunay tiling. In the case of the 5-fold symmetric 
quasicrystal, the optimum has wall width w/a = 0.103; 
the TE band gap is Auj/ujc = 42.3%, the largest ever 
reported for a photonic quasicrystal. For the octagonal 
quasicrystal, the optimal structure has dielectric mate- 
rial width w/a = 0.106 along each edge, which produces 
a TE band gap of Auj/uj c = 39.2%. 

The TE band gap formation is analogous to the TM 
case in that it involves an interplay between scattering 
from individual cells and the Bragg scattering of the 
quasiperiodic arrangement of scattering planes. Similar 
to the case of TM polarized radiation, we also obtain 
that the band gap for the optimal structures always oc- 
curs between bands Np and Np + 1 where Np is now the 
number of cells in the structure. 
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FIG. 5: Panels (a)-(d) show the electric field distribution for 
a 5-fold symmetric photonic quasicrystal optimized for TM 
polarized radiation. The structure consists of dielectric cylin- 
ders of radius r/a = 0.177 placed at the vertices of a Penrose 
tiling and displays a TM photonic band gap of 36.5%. Lower 
(a) and upper (c) band edge modes display a well-defined de- 
gree of localization; modes just below the lower band edge (b) 
and just above the upper band edge (d) display an extended 
character. Panels (e)-(h) show the magnetic field distribution 
in 5-fold symmetric quasicrystalline network optimized via the 
protocol in Sec. |Il]for TE polarized radiation. The structure 
consists of trihedral network with wall thickness w/a = 0.102 
and displays a TE PBG of 42.5%. The lower (e) and upper 
(g) band edge modes display a high degree of spatial concen- 
tration and modes just below the lower band edge (f ) and just 
above the upper band edge (h) display an extended character. 



FIG. 6: Panels (a)-(d) show the electric field distribution in 
8-fold symmetric structures for TM polarized radiation. The 
structure consists of dielectric cylinders of radius r/a — 0.189 
and displays a TM PBG of 42%. Lower (a) and upper (c) 
band edge modes displaying a well defined degree of localiza- 
tion; modes just below the lower band edge (a) and just above 
the upper band edge (d) displaying an extended character. 
Panels (e)-(h) show analogous magnetic field distributions in 
8-fold symmetric heterostructures optimized via the protocol 
in Sec. Ultbr TE polarized radiation. The structure consists 
of trihedral network with wall thickness w/a = 0.101 and dis- 
plays a TE PBG of 40%. Lower (e) and upper (g) band edge 
modes display a well-defined degree of spatial concentration; 
modes just below the lower band edge (f) and just above the 
upper band edge (h) display an extended character. 



Complete Band Gap: For the complete band gap, we 
find that the optimal structure consists of placing a wall 
along each edge and a cylinder at each trihedral vertex 
of the network generated by the protocol in Sec. HB The 
optimum for the 5-fold symmetric case has cylinder ra- 
dius is r/a = 0.157 and cell wall thickness w/a = 0.042. 
For the 8-fold symmetric point patterns, the optimum 
is r/a = 0.167 and w/a = 0.014. The scattering prop- 
erties of the individual scattering centers and cells are 
again essential in the band gap opening and the com- 
plete band gap occurs always between bands 3Np and 
3Np + 1, where Np is the number of points in the pe- 
riodic approximant point pattern (there are 3Np total 
scattering units in the system, 2Np dielectric disks and 
Np dielectric cells). 

The resulting optimal 5-fold symmetric structure dis- 
plays complete (TM and TE) PBG of 16.5% — the first 
complete band gap ever reported for a photonic qua- 
sicrystal with 5-fold symmetry and comparable to the 
largest band gap (20%) found for photonic crystals with 
the same dielectric contrast. 27 The optimal 8-fold sym- 
metric structure has a full PBG of 13.5%. 

As shown in Fig. 21 the band gap size obtained by con- 
sidering a contour along the Brillouin zone is completely 
consistent with the density of states calculations. Simi- 
lar to the case of the TM and TE PBGs, the system-size 
study presented in Fig. [3] reveals that the band gap size 
becomes independent of the size of the periodic approx- 
imant when Np > 100 (with Np the number of points 
in the original point pattern), implying that our results 
have converged to the large-system limit. 

For a photonic crystal, the states at the photonic band 
edges are propagating modes such that the electromag- 
netic field is distributed throughout the system. If the 
periodicity is disturbed, localized states begin to fill in 
the gap so that the states just below and just above be- 
come localized. Although the nature of the electromag- 
netic states in quasiperiodic structures needs to be fur- 
ther investigated, similar features are displayed by pho- 
tonic quasicrystalline systems. In Figs. [5] (a)- (d) and [6] 
(a)-(d), we compare the azimuthal electric field distribu- 
tion for modes below or above the band gap (plots (b) and 
(d) in Figs. |5] and |BJ), which display an extended charac- 
ter with the field distributed among many sites; and then 
modes at the band edges (plots (a) and (e) in Figs. [5] and 
[6]), which show a high degree of spatial concentration. 
We also find that the formation of the TM band gaps is 
closely related to the formation of electromagnetic reso- 
nances localized within the dielectric cylinders (as shown 
by the (a) and (b) plots in Figs. [5] and [6)) and that there 
is a strong correlation between the scattering properties 
of the individual scatterers (dielectric cylinders) and the 
band gap location. In particular, the largest TM gap 
occurs when the frequency of the first Mie resonance co- 
incides with the lower edge of the photonic band gap J£ 
Analogous to the case of periodic systems, the electric 
field for the lower band-edge states is well localized in the 
cylinders (the high dielectric component), thereby lower- 



ing their frequencies; and the electric field for the upper 
band-edge states are localized in the air fraction, increas- 
ing their frequencies. As shown in Figs. [5] (e)-(h) and [6] 
(e)-(h), an analogous behavior occurs for the azimuthal 
magnetic field distribution for TE modes: for states near 
the lower edge of the gap, the azimuthal magnetic field 
is mostly localized inside the air fraction and presents 
nodal planes that pass through the high index of refrac- 
tion fraction of the structure, while states near the upper 
edge display the opposite behavior. 



IV. CONCLUSIONS 

In sum, we have shown that it is possible to expand 
the spectrum of dielectric materials with sizeable com- 
plete band gaps to include photonic quasicrystals by in- 
troducing a novel constrained optimization procedure. 
The quasicrystalline photonic structures can be manu- 
factured using standard fabrication techniques used for 
photonic crystals^ 

Photonic crystals, the quasicrystals considered here, 
and the disordered heterostructures discussed in Ref. [l8| 
are all decorations of hyper uniform patterns^ A point 
pattern is hyperuniform if the number variance within 
a spherical sampling window of radius R (in d dimen- 
sions) grows more slowly than the window volume for 
large R, i.e., (N%) - (N R ) 2 = AR P , where p < d. Two 
dimensional crystal and quasicrystal patterns both cor- 
respond to p = 1. Our conjecture in Ref. [l8|, which was 
supported by a variety of examples, stated that a higher 
degree of hyperuniformity (smaller coefficient A) is ad- 
vantageous for obtaining substantial complete PBGs. In 
the Appendix, we provide the first calculations of the 
coefficient A for two-dimensional quasicrystals, namely, 
the 5- and 8-fold examples considered here. The num- 
ber variance for quasicrystals is greater than for crys- 
tals (for the same density) and less than the variance of 
disordered heterostructures. The results here are con- 
sistent with our conjecture in that the width of the op- 
timal complete band gaps for the photonic quasicrystals 
lies between the optima found for crystals and disordered 
heterostructures . 

Although photonic crystals have slightly larger com- 
plete band gaps, quasicrystalline PBG materials offer 
advantages for many applications. In the case of qua- 
sicrystalline structures the PBGs are significantly more 
isotropic, which is advantageous for use as highly-efficient 
isotropic thermal radiation sources^ The properties of 
defects and channels useful for controlling the flow of light 
are different for crystal and quasicrystal structures. Ra- 
diation with frequencies above or below the band edges 
are propagating modes that are transmitted through pho- 
tonic crystals but are more localized modes in the case 
of quasiperiodic patterns, which give the former advan- 
tages in some applications, such as light sources^ On the 
other hand, due to the lack of translational symmetry, 
PBG structures constructed around quasiperiodic point 



patterns can provide a large number of inequivalent local 
environments and as such can support a rich variety of 
localized modes. These localized modes can be employed 
to design laser systems with highly unusual field patterns 
with possible applications to biological sensing. 31 
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Appendix: Asymptotic Number Variance for the 5- 
and 8-Fold Quasicrystal Point Patterns 

Here we compute the asymptotic number variance for 
the 5- and 8-fold quasicrystal point patterns in two 
dimensions obtained from the Penrose and octagonal 
tilings. For such point patterns, the number variance for 
large windows grows like the perimeter of the window, 
i.e., 



(N 2 R )-(N R f=A^ + 0(l), 



R 



oo, 



(A.l) 



where R is the window radius and D is the mean-nearest 
neighbor distance associated with the point pattern. In 
order to compare the dimensionless surface area coeffi- 
cient A for different point patterns, the comparison must 
be made at the same number density p (number of points 
per unit volume). It was shown in Ref. 28 that this 
amounts to rescaling A such that A/y/cj), where 



7T ^o 



(A.2) 



is an effective packing fraction. Thus, we define the 
rescaled coefficient A = A/y/cj). 

The rescaled surface area coefficients can be evaluated 
numerically by employing finite large quasicrystalline 
point patterns (of around 15,000 points) and count the 
number of points in a circular window of radius R. For 
a fixed radius, the circular window scans the quasicrys- 
talline domain and the variance of the number of points 
in the window is evaluated (the windows employed in 
our calculations contain from a few to 5000 points and 
to calculate the variance we use 1000 configurations with 
a fixed window size). The surface area coefficient is then 
obtained by analyzing the behavior of the number vari- 
ance in the limit of large window size. 

We find that the rescaled surface area coefficients for 
the 5- and 8-fold case are given respectively by As.foid = 
0.60052 and Ag-foid = 0.59567. Our results for these qua- 
sicrystal cases should be compared to results for some 
periodic point patterns; for the triangular, square, hon- 
eycomb and Kagome lattices, A tr i an guiar = 0.508347, 

^-square = 0.51540, ^-honeycomb = 0.55702 and ^Kagome = 

0.58699, respectively^ It is noteworthy that the triangu- 
lar lattice has been proven to have the smallest surface 
area coefficient among all lattices and is thought to be 
the global minimum among all point patterns. Thus, 
we see that the quasicrystal point patterns studied here 
have a high degree of hyperuniformity, even if not as large 
as that for the aforementioned crystals. These calcula- 
tions of the surface area coefficients and PBGs reported 
here are consistent with the conjecture made in Ref. 
1181 . namely, that the width of complete PBGs of high- 
dielectric contrast photonic structures is correlated with 
the degree of hyperuniformity of the underlying point 
pattern. 
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